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Abstract 

We solve the complex extension of the chiral Gaussian Symplectic Ensemble, defined as a Gaus- 
sian two-matrix model of chiral non-Hermitian quaternion real matrices. This leads to the appear- 
ance of Laguerre polynomials in the complex plane and we prove their orthogonality. Alternatively, 
a complex eigenvalue representation of this ensemble is given for general weight functions. All 
fc-point correlation functions of complex eigenvalues are given in terms of the corresponding skew 
orthogonal polynomials in the complex plane for finite-A'^, where N is the matrix size or number 
of eigenvalues, respectively. We also allow for an arbitrary number of complex conjugate pairs of 
characteristic polynomials in the weight function, corresponding to massive quark flavours in appli- 
cations to field theory. Explicit expressions are given in the large- limit at both weak and strong 
non-Hermiticity for the weight of the Gaussian two-matrix model. This model can be mapped to 
the complex Dirac operator spectrum with non- vanishing chemical potential. It belongs to the sym- 
metry class of either the adjoint representation or two colours in the fundamental representation 
using staggered lattice fermions. 



1 Introduction 



The first complex random matrix models date back to Ginibre ^ where the three classical Wigner- 
Dyson ensembles were generalised. In recent years we have seen a revival of complex matrix models, 
due to new applications as well as new insights and solutions. Todays interest in such models ranges 
from Quantum Chromodynamics (QCD) with chemical potential j2|, fractional Quantum-Hall effect 

Coulomb plasma and growth processes |E] to scattering in Quantum chaos reviewed in [Hj. 
New insights include the discovery of the regime of weak non-Hermiticity [7], interpolating between 
correlations of real eigenvalues of the Wigner-Dyson ensembles and those correlations computed by 
Ginibre. Furthermore, we know that there are many more complex matrix model symmetry classes 
than those with real eigenvalues jHl. 

The aim of this work is to solve such a new complex matrix model, by generalising a chiral version 
of one of the Wigner-Dyson ensembles. Although this is motivated by a specific application to QCD 
we hope that our result will be useful in other areas as well, because of the wide range of applicability 
of matrix models. A major problem in QCD is that the Dirac operator becomes non-Hermitian when 
introducing a chemical potential for the quarks. In general this spoils the positivity of the action and 
thus its interpretation as a Boltzmann weight when performing numerical simulations of QCD on an 
Euclidean space-time lattice. 

A way out of this dilemma is to study properties of field theories similar to QCD, where Dirac 
eigenvalues come in complex conjugate pairs. There are two symmetry classes different from QCD, 
with the Dirac operator having an orthogonal or symplectic symmetry compared to unitary. This 
feature holds at both zero [Hj and non-zero chemical potential JHI- In this work we will construct 
a complex matrix model with symplectic symmetry and solve it analytically. It has the virtue to 
be testable in lattice simulations including dynamical fermions. These are absent in the so-called 
quenched approximation, suppressing the Dirac operator determinants, and such matrix model pre- 
dictions im I12j have already been confirmed in quenched lattice simulations of QCD and the 
symplectic symmetry class |14| . where preliminary results of this work were used. Let us emphasise 
that unquenched matrix model predictions for QCD with chemical potential do exist jl5l I16j . 

The most important ingredient of matrix models is the way they incorporate the repulsion between 
eigenvalues or levels. The form of the repulsion term is dictated by the Jacobian when changing from 
matrices to eigenvalues. For real eigenvalues it is given by the absolute value of a Vandermonde 
determinant of the eigenvalues to the power of the Dyson index (3 = 1,2, or 4, labelling the number of 
independent matrix elements in the Orthogonal, Unitary or Symplectic Ensembles respectively (OE, 
UE or SE). The same is true for the squared singular values of the so-called chiral ensembles (chOE, 
chUE or chSE) which have an additional block structure. When employing the method of orthogonal 
polynomials all eigenvalue correlations of the corresponding Gaussian ensembles (CUE etc.) can 
be obtained using Hermite polynomials, or Laguerre polynomials in the chiral classes. 

What is known about the level repulsion of complex eigenvalues and the corresponding correla- 
tions functions? For the complex CUE one obtains again the modulus squared of the Vandermonde 
determinant of complex eigenvalues For the corresponding chiral ensemble it was conjectured 
in ^2 to have the same Jacobian of squared complex eigenvalues, which was later obtained from a 
two-matrix model realisation ^2] ^ • While for /3 = 4 a Vandermonde to the power 4 can be constructed 
using normal matrices jl8l E| j the Jacobian of the complex GSE of generic non-Hermitian quaternion 
real matrices is more complicated It shows an additional level repulsion from the real line as 
pointed out in (see also [10]), where this ensemble was solved for an arbitrary weight function. A 
similar feature was observed numerically in [lOj . where the level repulsion for all three complex chiral 
ensembles were compared. 

^Chiral complex models with a single matrix were first defined in I^IIU] without having an eigenvalue representation. 
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Our task here is to derive the corresponding Jacobian for the complex chGSE using a two-matrix 
model of non-Hermitian quaternion real matrices, in analogy to the approach ^Sj. Because of the 
QCD application we allow for exact zero-eigenvalues, using rectangular matrices. The result for the 
Jacobian turns out to be again that of the non-chiral complex GSE of squared variables (generalised 
to rectangular matrices), hence having an additional level repulsion from the imaginary axis. 

We solve our model introducing orthogonal polynomials in the complex plane, where we encounter 
complex Laguerre instead of Hermite polynomials ^3,, 21, 19 in non-chiral ensembles. Other techniques 
include supersymmetry, e.g. in |22[ I2H I2()j . the replica method |23| I12j or a fermionic mapping 
|24j which we do not discuss here. Our computation completes the task of solving the complex 
extensions of the chiral and non-chiral Wigner-Dyson (3 = 2 and 4 ensembles. The difficulty with the 
remaining (3 = 1 ensembles is that they have a finite fraction of real eigenvalues, with the corresponding 
probabilities to be determined. Only the spectral density at weak non-Hermiticity of the complex GOE 
is know so far |22j . 

The article is organised as follows. In the following section [21 the complex chSE is defined for a 
general weight function, and all correlation functions are given at finite- in terms of the corresponding 
skew orthogonal polynomials. The subsections 12.11 and 12.21 distinguish between weights without and 
with explicit insertions of characteristic polynomials or mass terms respectively, where we have to 
impose a two-fold mass degeneracy. In subsection 12.31 two examples for weight functions and their 
corresponding skew-orthogonal polynomials are given, a Ginibre type weight and a Gaussian combined 
with a X-Bessel function. The latter originates from the Gaussian two-matrix representation of our 
model discussed in detail in section [3 The orthogonal polynomials of the ii'-Bessel function weight 
are Laguerre polynomials in the complex plane, and we prove both their orthogonality and skew 
orthogonality, with the help of appendix EF- In section |2l we take the large- limit at weak non- 
Hermiticity for the ensemble of complex skew orthogonal Laguerre polynomials. Explicit examples 
for the quenched and massive spectral density as well as for massive partition functions are given in 
subsections 13.11 and 13.21 . Section |1] is devoted to the strong non-Hermiticity limit, giving the same 
set of examples. Technical details for the derivation of the limiting kernel at strong non-Hermiticity 
are collected in the appendix IbI In section |3 we relate the Gaussian two-matrix model representation 
of our complex model to the Dirac spectrum with chemical potential. We also provide a matrix 
representation of characteristic polynomials. The derivation of the Jacobian is deferred to appendix 
lO Our findings are summarised in section |21 

2 Correlation functions of the complex chSE at finite- 

In this section we define the complex extension of the chSE for a general class of weight functions. 
The diagonalisation of the matrix model representation in terms of two non-Hermitian quaternion real 
matrices is discussed in more detail in section El and appendix O This Gaussian matrix model leads 
to a specific weight function containing i^-Bessel functions. The reason is the nontrivial decoupling of 
eigenvalues and eigenvectors which is not present for ensembles with real eigenvalues. The two-matrix 
model representation determines the Jacobian we write down below as the interaction term for the 
eigenvalues. 

For the calculation of the eigenvalue correlation functions in terms of skew orthogonal polynomials 
we proceed in two steps. We first treat the so-called quenched case where no explicit characteristic 
polynomials (determinants) are present in the weight function. Although the expressions we give 

■^The Gaussian weight first introduced in for complex Laguerre polynomials LJv agrees only asymptotically or 
at special values of v with the correct _Rr-Bessel weight. The latter was introduced in ^ 15 and conjectured to satisfy 
orthogonality. 
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for the complex eigenvalue correlations in the first step, eq. (|2.7j) below, are true for general weight 
functions, it is easier to include determinants as mass terms in a second step. Namely we can re- 
express the corresponding massive correlation functions in terms of the quenched ones, following the 
same logic as in |^ . 



2.1 Complex eigenvalue correlation functions 



The matrix model partition function of the complex Gaussian chSE is defined as the following matrix 
integral 



d^d^ exp 



-iVTr($t$ + xft^j 



(2.1) 



where $ and ^ are two matrices with quaternion real elements of rectangular size (A^ + z^) x without 
further symmetry properties. In order to transform to eigenvalues we parametrise the following linear 
combinations as 



C = = U{X + R)V 

D = i^t + ^^t = v\Y + S)U^ 



(2.2) 



being equivalent to independent Schur decompositions oi CD and DC. U and V are symplectic 
matrices and the diagonal matrices X and Y contain the complex conjugate eigenvalue pairs {xk,x*^) 
and {yk,yk), respectively. After integrating out the second set we obtain in terms of the variables 
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the following expression for the partition function 



N 



N 



2|2 



4 



* 2|2 



k>l 



N 

n 

h=l 



^*2|2 



(2.3) 



(2.4) 



For more details we refer to section [S] and appendixinjfor the Jacobian. The integrals in Zj. € C extend 
over the whole complex plane with d^z = d^e{z)d^m.{z). Here w{z.,z*) denotes a general positive 
definite weight function defined in complex plane. It depends on both z and z* . The only restriction 
we want to impose here is the existence of all positive moments, J d'^z w{z, z*) z'^^ < oo, k £ N , and 
similarly for the complex conjugated ones. For the Gaussian matrix model eq. (|2.1jl the weight is 
given by (see section 13)) 



w 



(2u) 
K 



iV(l + ^^ 



2^^ 



exp 



/ie (0,1] 



(2.5) 



a second example is given below in eq. (|2.27|) (see also [27,). The general weight w{z, z*) may depend 
on extra parameters such as the non-Hermiticity parameter /i > 0, or masses mj as we will discuss 
later. The former will allow us to take limits to either real eigenvalues by letting ^ ^ 0, or to maximal 
non-Hermiticity at — > 1, given for example by the chiral Ginibre type weight |z| exp[— zz*]. Eq. 
()2.4() can be derived in a different way, inserting squared arguments into the Jacobian of the non-chiral, 
complex SE This is in analogy to the construction of the chiral from the non-chiral complex UE 

Below we compute all eigenvalue correlations in terms of skew orthogonal polynomials in the 
complex plane. This generalises the construction using complex skew orthogonal polynomials for 
the non-chiral complex SE. There, complex Hermite polynomials were used for the special case of a 
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Gaussian potential. In our chiral model with the special weight of the Gaussian two-matrix model we 
encounter complex Laguerre polynomials, as expected from the complex chUE (see 115| ). 

For a general weight function the A;-point correlation functions of complex eigenvalues can be 
defined in the standard way 



Rn{zi, ...,Zk) 



. . N N N N 



(N-k) 



^ j=k+l i=l n>l h=l 



Being real positive functions they depend also on the complex conjugated variables zj, . . . , z^. We 
suppress this dependence in the following for simplicity. In analogy to the non-chiral model the 
eigenvalue correlations show a depletion of the eigenvalues from the real axis, being identically zero 
there. This property was found also numerically for a chiral complex one-matrix model with chemical 
potential ^Hl- Since our Jacobian or interaction term is given in terms of squared variables compared 
to jn]) there is the same depletion also from the imaginary axis, due to — z)^^ = 4:i^e{zh)Qm.{zh)- 
Such a depletion has been already detected in lattice data for the corresponding symmetric class [Tl] . 

Following 1171 119j all fc-point correlation functions are determined through a quaternion determi- 
nant 17^ of a matrix valued kernel, 

RN{zi,...,Zk) = Qde\j^i^_^f,[)CN{zi,Zj)] , (2.7) 

where the kernel is given by 

JCn{zi,Z2) = -^1^2 - Z2)w{zi,Zi)w{Z2,Z2)\' [ ^ . ^ , Jn • (2-8 

V I'^N{Zi,Z2) -Kn[Zi,Z2) J 

Here we have introduced the pre-kernel kn{zi,Z2), ^ 

KNizi,Z2) = ^2 — iq2k+iizi)q2kiz2) - q2k+iiz2)q2k{zi)) ■ (2.9) 
k=Q 

Examples for such correlation functions are the spectral density given by 

Rn{z) = [z*"^ - z^) w{z,z*) i^n{z,z*) , (2.10) 

and the two-point function 

Rn{z,u) = {z*'^ — z'^){u*'^ — u^) w{z, z*)w{u,u*) 

X [KN{z,Z*)KMiu,U*) -\kn{z,u)\'^ + \KNiz,U*)\'^) . (2-11) 

The pre-kernel eq. (|2.9j) is expressed in terms of skew-orthogonal polynomials qk{z)- These are 
polynomials of degree k in the complex variable z G C squared, qk{z) = z"^^ + 0{z'^^^'^), where we 
have chosen the monic normalisation (because of the form of eq. (|2.4|) we only need polynomials of 
squared arguments, see also ITT]). The skew-orthogonal polynomials are defined with respect to the 
following antisymmetric scalar product 

{h\g)s ^ [d'zw{z,z*){z*^-z'){h{z)g{zr-h{zrg{z)), (2.12) 



^For simplicity we do not give the more general expression of )19| here. 
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by obeying 



{<l2k+l\<l2l) S = -{Q.2l\<l2k+l)s = rk Ski , 

{q2k+i\q2i+i)s = {q2i\q2k)s = o . (2.13) 

For a general weight they enjoy the following representation 

q2N{z) = /ll{z'-zf){z'-z;')\ , (2.14) 

IN N \ 

q2N^i{z) = IT\{z'-z^){z^-z*,'){z^ + Y,{4 + 4') + cn)) , (2.15) 

\l=l k=l I Zm 

which holds in complete analogy to T^. Here cat is an arbitrary constant that can be set to zero. Due 
to this representation we can observe that the polynomials have real coefficients, as q2N{z)* = q2N{z*) 
and q2N+i{z)* = q2N+i{z*). A similar representation of the skew orthogonal polynomials also holds 
for the SE with real eigenvalues |28j . In section [S] we will also give a matrix representation of the 
q2N{z) in terms of a characteristic polynomial. The expectation values are to be taken with respect to 
the partition function eq. ()2.4() (not to be confused with the antisymmetric scalar product). The skew 
orthogonal polynomials q2N[z) and q2N+i{z) both depend only on the variable z. While eq. (|2.7() for 
the A;-point correlation functions holds for general weight functions it is convenient to explicitly split 
off determinants or mass terms from the weight. We shall do this in the next subsection, allowing us 
to compute characteristic polynomials more general than in the example eq. (j2.14|) above. 



2.2 Massive complex correlation functions 

The massive matrix model partition function of the complex chSE we wish to compute is given by 

„ iV Nf 
Z^^'^'himj}) = n{Sz,w{z,,z*)\{\rnfr\z^^+m)\\^^+mf\'' 

N N 

x\{\4-zl?\zl-zt'\^\{\zl-zl^\\ (2.16) 

k>l h=l 

where we allow for complex masses 7^ m^.^ is this subsection. Compared to the partition function 
eq. 1)2. 4|) we have explicitly introduced mass terms in complex conjugated pairs as extra sources. The 
new weight function including the masses is still positive definite^. The mass terms can also be written 
as a determinant of quaternion real non-Hermitian matrices as we will show in section |21 Therefore 
these partition functions can be interpreted as characteristic polynomials, and we will give explicit 
formulas for them below. 

In the presence of a non-Hermiticity parameter ^ we can make contact with the chSE by taking the 
Hermitian limit ^ 0. This leads to 4-fold degenerate (real) masses, as indicated by the superscript^. 
Our aim is to express the partition function as well as all complex eigenvalue correlation functions in 
the presence of mass terms through the eigenvalue correlations without mass terms at Nf = 0. For 
that reason we need to impose such a fourfold mass degeneracy as in [2^. There, massive correlation 

''This is no longer true for mass terms in the complex chUE corresponding the the QCD symmetry class. 
*Due to Kramer's degeneracy the mass terms automatically come in doubly degenerate pairs for the symplectic 
ensemble, see also sect. |K1 
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functions were expressed through quenched ones for Hermitian self-dual matrices [22], and we will 
follow the same strategy here. Similar results were also obtained for complex non-Hermitian matrices 
in the non-chiral case (3 = 2 jSHj, where only a twofold degeneracy is needed due to the different 
Jacobian H^^^ \zl-zf\^. 

To start we compute the massive partition function eq. ()2.16() in terms of quenched correlation 
functions. When writing the mass terms as z? + = -^j — (imj)^ they can be absorbed into a larger 
Jacobian of + Nf eigenvalues zi,. . . ,zi\f,imi, . . . ,imj\fj,, after multiplying with constant factors as 



m 



* 2|2 
/ 



7^ 0. The massive partition function is thus proportional to a quenched A'^^-point 
correlation function at i times the complex masses. Filling in all explicit factors from the definition 
eq. (|2.6() we thus obtain 



-'N 



{{mj}) 



m 



{N + Nf)] Yl^Li w{imf, —im*j:)\m'j — m*^ 



* 2|2 



m 



2|2 



mt 



m 



* 2|2 



Nf 



Au 



■'N+Nf 11/: 



Qdetj^k=i,.. .,Nf 



KM{-im*,imk) 
K]y{imj,imi^) 



-KN{-im*, -iml) 



Nf 
k>l 



m 



212 



mt 



m 



* 2|2 



(2.17) 



Dividing by Z^^_^^^ and switching back to parameters uj = imj this gives an explicit expression for 
products of L pairs of characteristic polynomials times their complex conjugates 

I N L \ 



\j=i 1=1 



* 2 



n^){z^ - ui'){z^ 



* 2 



Av 



Qdet 



Lk=l, 



KN[Ui,Uk) 



' Zn 

-kn{u*,uI) 
-kn{ui,uI) 



{N + Ly. Yl 



n 



k>l I""! 



U 



2|2| 



U 



* 2|2 



(2.18) 



It extends the corresponding results for the corresponding (3 = 2 class. Next we use the same 
trick to compute the massive eigenvalue correlation functions, obtaining 



,Zk) 



{N + Nf)l 4^^/)({^^}) 



(2.19) 



R 



■j^+Ar^(^i, • • .,Zk,imi, . . .,imNj] 



n^/ I- 



1/ 



-imj)\mj 



m 



*2|2 



Nf I 2 

hii IK 



m 



212 



mt 



m 



* 212 



Inserting the above relation ()2.17|) for the massive partition function all pre-factors cancel and we 
obtain the very elegant expression 



R 



(ANf) 



N 



{Zl, ...,Zk) 



R 



R 



(0) /. . N 



(2.20) 



7 



The same relation was obtained previously for real eigenvalues in the chSE and for the non-chiral 
complex UE j^H]- It is easy to see that it also holds for the non-chiral complex SE, replacing squared 
by un-squared variables everywhere. Together with eqs. p.7j) and (|2.8|) this gives all eigenvalue 
correlation functions for our massive model eq. ()2.16|) . Due to the different numbers of eigenvalues 
occurring, A'^ compared to N+Nf, these relations are only exact for an A^- independent weight function 
w{z, z*). When we consider A^-dependent weights as in the examples eqs. (|2.27|) and 1)2. 21|) below our 
results will be still valid for asymptotically large- as then N ^ N + Nj in the exponent. 

We have to add a word of caution to applying eqs. (|2.2()j) and (|2.17)) . While these expressions are 
true for general complex masses ruf the limit of real masses has to be taken with care. As we have 
explained after eq. (|2.6|) our model shows a depletion of eigenvalues both from the real and imaginary 
axis. Therefore both the numerator and the denominator in eq. ()2.2U|) vanish when setting the masses 
TUf to be real or purely imaginary. We thus have to take the limit of the imaginary part of going 
to zero using the rule of de I'Hospital. This leads to a finite expression for the massive correlation 
functions. We will illustrate this explicitly in examples in sect. |31and0when taking the large- A^ limit. 



2.3 Examples for skew orthogonal polynomials in the complex plane 



In order to illustrate the finite-A^ expressions for the eigenvalue correlation functions we give two 
explicit examples of weight functions together with their corresponding skew orthogonal polynomials 
eqs. 1)2. 14() and 1)2. 15(1 as they appear inside the pre-kernel. We will use some results for non-skew 
orthogonal polynomials with respect to the same weight functions which are derived in appendix 
The first example is from the Gaussian two-matrix model representation of the complex chSE (see 
sect. ISJ, 



w 



i2u) 
K 



iz,z*) 



A^(l + /i^)|z|' 

2^2 



exp 



{z^ + z 



* 2-y 



;UG(0,1]. (2.21) 



Here is a non-Hermiticity parameter. Compared to the QCD symmetry class |15| the parameter v 
measuring the rect angularity of the matrices is shifted, v — > 2z/, as for the real eigenvalue model. The 
real limit can be taken by letting /u ^ 0, leading to a delta- function in the imaginary part. This will 
be discussed in more detail in sect. 01 below. We only mention here that in this Hermitian limit the 
last term of the Jacobian in eq. 1)2. 4() will provide two more powers of the real part, (5Re(2;))2, arriving 
at the standard power of zero-modes |5Re(2;)|^'^+^ for the real chSE (see eq. (|5.1flj) ). 

Let us consider a special case in which the ET-Bessel function simplifies. For indices v 



exactly coincides j29| with its asymptotic value, K2u{x) ~ e 



w 



(2u=±\) 
K 



{Z,Z* 



1 



exp 



-A^(l + M^ 

2/i2 



+ 2- 



2x ' 



Ifl 



±i it 

and we obtain up to a constant 



1 + ^2 



iz^ + z*^ 



(2.22) 



This is of the same form as the weight in the corresponding non-chiral ensemble |19| . More generally 
speaking the i^T-Bessel function at half-integer index simplifies to an exponential times a finite sum over 
half-integer powers of its argument. In the non-chiral complex SE _19_ complex Hermite polynomials 
were used to construct the skew orthogonal polynomials eq. 1)2. 13|) for a Gaussian weight function as 
eq. p.22|) . This is consistent with the relation between odd or even Hermite polynomials and Laguerre 

polynomials , respectively (see e.g. [221) • 

Following the same lines as in jl9j we construct our skew orthogonal polynomials here using 
complex Laguerre polynomials, which are show to be orthogonal with respect to the weight eq. 1)2.21(1 
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in appendix EI We obtain 



qg^,{z) = -{2k + l)l(^ljfY^\ll^,(^^y (2.23) 

q,,{z} - /C. i (fc + ZV + lj Z.(l + ^2)2,22j-j-!r(i + l.+ l)^2, Vl_^2 

(2.24) 

in monic normalisation q2k+i{z) = (z^)^^^-'^ + . . . and q2k{z) = {z"^)"^^ + . . . 

The fact that this choice satisfies eq. (j2.13|) can be seen as follows. In the case where both 
indices in eq. ()2.13|) are equal, the second line for k = I, the scalar product vanishes trivially due 
to antisymmetry. If both indices are odd we now show that the scalar product vanishes because of 
the orthogonality of the Laguerre polynomials from appendix To see that we use the three-step 
recursion relation 

Z'Ll^iZ') = -{n + l)Ll\,{Z^) + {2n + 2v + l)Ll%Z^) - (n + 2u) LI^_,{Z') , (2.25) 

in the squared variable Z"^ = jz^- Consequently the factor (z*^ ~ z"^) the scalar product eq. 
(|2.12|) can only lower or raise the index of the Laguerre polynomials by one unit. The terms where 
indices are raised or lowered vanish from orthogonality, having different parity. The middle term in 
the recursion eq. (|2.25p only contributes at k = I and the corresponding terms cancel (as we have 
seen already from antisymmetry above). The same logic applies if both indices in eq. (|2.13j) are even. 
Only those terms in the recursion where the index is unchanged contribute, and they also cancel after 
using orthogonality, for all values of k and /. 

It remains to analyse the scalar product in eq. ()2.12|) with one index even and one odd, {q2k+i,Q2i)s- 



For k > I the expression again vanishes trivially due to orthogonality, using eq. (|2.25j) . So far the 
coefficients in the linear combination of even Laguerre polynomials in the definition 1)2.24(1 could have 
been arbitrary. It is the requirement {q2k+i, q2i)s = fk^ki for k < I that fixes them uniquely in a given 
normalisation. One can easily check that the coefficients given in the definition (|2.24|) imply vanishing 
for k < I, upon using the recursion and the norms of the Laguerre polynomials eq. ()A.19|) . The overall 
constants in eqs. ()2.23() and (|2.24|) are fixed by choosing the monic normalisation. This choice allows 
us to take the limit ^ 1 of maximal non-Hermiticity on the polynomials. The remaining coefficients 
in eq. ()2.13|) at A; = / follow to be 

rf = 8V(2fc + 1)! r{2k + 2u + 2)^-^£^L__ . (2.26) 

This determines all quantities in the pre-kernel eq. (|2.9p and thus all A;-point correlation functions 
through eq. (|2.7|) . 

We now turn to the second example given by a Gaussian, chiral Ginibre type weight 

w^^^^z^z*) = Izl^'^+iexp [-7V|z|2] . (2.27) 



In can be obtained from the above weight eq. ()2.21|) at maximal non-Hermiticity ^ = 1, when taking 
the asymptotic limit of large arguments (or at v = ±5)- However, at small arguments they differ, 
and so far no matrix representation is known for the weight eq. (|2.27() . As for all rotational invariant 
measures the associated orthogonal polynomials are the monomials z^ and their complex conjugates, 
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from which we can construct the skew orthogonal polynomials. We obtain the following result 



T{k + - + - 



1 +2 + 8;^24.r(, + ^ + |)r(j + | + |) ' 



(2.28) 



(2.29) 



in monic normalisation. The skew orthogonality eq. (|2.13p can be checked along the same lines as 
described above for the Laguerre polynomials, with a much simpler recursion of course. We only give 
the final result for the coefficients r^, using the normalisation of the orthogonal polynomials computed 
in the appendix lAleq. (|A.22|1 . 



r(4A; + 2i/ + I) 



27r— ^ . (2.30) 

Let us also give an example for the simplest correlation function at finite- A^, the spectral density. 
Taking the simplest case without determinants and the weight eq. (|2.21|1 of the two-matrix model we 
can explicitly write out eq. (|2.10|) : 



Rn{z) 



2 



^2\u|4i^+2 7^ 



iV(l + ^^)|2| 

2^^ 



exp 



{z^ + z 



* 2^, 



(2.31) 



iV 



8^/^"^ ^ (1 + /x2)2fc+2j+2i/ p(2fe 2i/ 2)j! r(j ^v^\) 



(1 - ^2)2fc+2i+i 22'=-2i^! r(A; + + l)(2j)! 



X L 



-2v 
^2k+l 



Nz 



* 2 



1^ 



L 



2v 

2fc+l 



1^ 



L: 



* 2 



It is manifestly real and vanishes along the real and imaginary axis. 



3 The large-A/^ limit at weak non-Hermit icity 

The limit of weak non-Hermiticity was first introduced in j7j for unitary invariant complex non- 
Hermitian matrices. It is defined by simultaneously taking the Hermitian limit /i — > and — > oo 
such that the following combination 

lim 2Nii^ = o? (3.1) 

is kept constant. In this limit the macroscopic spectral density collapses onto the real axis, while the 
microscopic correlation functions still extend into the complex plane, on a stripe of the order of o? . 
This limit allows us to extrapolate between the correlation functions of real eigenvalues, by taking 
a ^ 0, and those of complex eigenvalues at strong non-Hermiticity to be introduced in the next section 
^ by taking a — > oo. We define the following rescaling of the complex eigenvalues at the origin 

\/2 iV(5Re z i9m z) = ^Nz =i , (3.2) 

where we take over the convention from /i = 0. There eigenvalues are rescaled as VYjX = ^ where 
2A^ = y is the volume, and the inverse variance or chiral condensate proportional to the macroscopic 
density at the origin is given by S = l/-\/2 in the matrix model eq. 1)5. 1(1 at = 0. The pre-kernel 
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and the resulting microscopic correlation functions have to be rescaled accordingly: 

1 ^ ?! ^2 



In the following we derive these quantities explicitly for the weight eq. (|2.2H) . 

We start with the large- weak non-Hermiticity limit of the skew orthogonal Laguerre polynomials 
eqs. 1)2. 23p and H2.24() . Here, we can use some of the results from which we repeat for completeness. 
From the asymptotic of the Laguerre polynomials, 

lim k-^^LfiZ^) = {kZ^)-''J2^ (2VkZ^) , (3.4) 

fc— >oo V / 

we obtain for the odd Laguerre polynomials of eq. 1)2. 23(1 

(^) = A'-2-,.'r-^.„ (2V-S 0.^-2,,. (3.5) 

Here we have introduced the scaling variable s G [0, 1]. For the even polynomials eq. (|2.24|) we have to 
replace the sum by an integral, X]j=o ~^ ^ lo introducing a second scaling variable t = j/k ^ [0, 1]. 
In terms of the integration variable t we obtain for the even Laguerre polynomial inside the integrand 

hm Lf (-I^) = N^-2'-{strC^^J2j2V^td , ^^t, (3.6) 
using 2jNz'^ — > tsS,'^. The normalisation constant is mapped to 

i,fc,iv^oo 22jj!r(j + + 1) ^ ^ ' 

The resulting expression for the sum over Laguerre polynomials eq. 1)2. 24p is then 

.,^J™.^o2-. [iT7) 22.,!r(, + . + i)^^^- I J - X ^l~t ' ^-(2V.tO. 

(3.8) 



j=0 



Here we have used the fact that the //-dependent pre-factors turn into an exponential, 

lim f(l±/.2)2^' = e2^''^(i±"'/2JV)\ ^ ^±sta^ (39) 



and similarly for powers of k. In order to obtain the pre-kernel eq. (|2.9I) it remains to evaluate the 
fc-dependent pre-factors of both skew orthogonal polynomials eq. (|2.24() and 1)2. 23p . divided by the 

k ' 



norm rj^ , in the large- limit. 



lini -l^(2k + iy :^~/2, k\T{k + u + l f ^\+/^ = le-2-^ (3.10) 

Having collected the asymptotic of all building blocks we can write down the limiting expression for 
the pre-kernel eq. (|2.9() 



1 9 



2u-3 



^'ds^';^e-2^(i+*)"' (j2.(2Vit6)J2.(2v^e2*) - J2u{2VTte2)J2u{2V~sCi] 
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The asymptotic of the even and odd skew orthogonal polynomials eqs. (|2.24|) and (|2.23j) follow from 
setting s = 1 in eqs. p.Hfl and ()3.5() times e^" , respectively. The former will be used later in 
subsection 13.21 as it enjoys a direct interpretation as a characteristic polynomial or massive partition 
function. 

In order to give the spectral density and all higher correlation functions we also need the limiting 
form of the weight eq. 1)2. 21j) : 



7V- 



lim 



W 



K ' 



Z, Z 



V 2a2 



exp 



1 

4a2 



* 2\ 



(3.12) 



Inserting this together with the limiting form of the pre-kernel eq. (|3.11jl into eqs. (|2.8jl and (|2.7jl 
rescaled according to eq. (|3.3() we obtain all quenched microscopic correlation functions, and through 
eq. 1)2. 20|) also all massive microscopic correlations functions. Since the massive correlation functions 
at real masses involve taking limits we give explicit examples in subsection 13.21 below, after first 
discussing the simplest example, the spectral density without masses. 



3.1 Example: the quenched spectral density 

Here we give the quenched [Nj = 0) microscopic spectral density for an arbitrary topological charge 
V as well as a comparison to its real limit a — > 0. It is given by inserting eqs. (|3.1H) and (|3.12|) into 
eq. (|TTn|l 



Pweak 



1 



32a4 



* 2 



e) K2. ) exp 



(3.13) 



Jo Jo yt ^ 

We note that in addition to the z^-dependence of the J-Bessel functions the weight is also explicitly 
topology dependent. This is in contrast to the spectral density of real eigenvalues, see eq. (|3.14j) 
below. It is only in the limit of small or large arguments, |i^p/(4a2) ^ 1, that the K-Bessel 



2x ' 



leading to a Gaussian weight asymptotically (see 



function becomes z^-independent, K2u\ 
also eq. ^TF^). 

In Fig. n we have plotted the microscopic density for = and 2 at the same value of weak 
non-Hermiticity parameter a = 0.4. It is instructive to look at the density of real eigenvalues x for 
the chSE ^ 



Preal{x) 



2x' 



duu^ / dv {J2v(2uvx)J2v+i{2ux) — vJ2v{2ux)J2v+i{2uvx)) 



(3.14) 



Compared to the real density the density in the complex plane keeps the oscillatory structure, with 
the maxima and minima at the same places. However, due to the repulsion of eigenvalues from the 
real axis the peak structure is doubled, with a valley along the real axis. On the right of Fig. I^lthe 
additional zero eigenvalues aX v = 2 repel the density further away from the origin, just as they do in 
the real case. The double peak structure is very different from the complex chUE |111 1121 115j and has 
been detected in lattice simulations jTl]. 

We can also perform an analytical check by taking the limit a — > on the microscopic density eq. 
(|3.13)1 . In this limit the pre- factor of the integral becomes 



lim 



1 



o32a4 



K2. [ ^ 



2ol- 



1 

64^ 



—x\Jx^ + y2 5(?/)' , 



16 



(3.15) 



12 



Pweak 



Pweak 




9m ^ 



Figure 1: The microscopic spectral density of complex eigenvalues at weak non-Hermiticity parameter 
a = 0.4 for 1/ = (left) and u = 2 (right). 
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Figure 2: The microscopic spectral density of real eigenvalues for v = Q (left) and v = 2 (right). 

where we denote ^ = x + iy. Taylor expanding the integrand in the density with respect to y = 9m(^), 
we obtain to leading order 

iim (j2,{2V7toj2u{2V^e) - (e^D) 

= Aiy^s(j2y{2^/7tx)J2u+i{2^sx)-Vij2u+i{2^/7tx)J2u{2^sx)] + 0{y^) , (3.16) 



after using a Bessel identity. For the limiting spectral density we thus obtain 

,7r 



a— >0 



4 

J ds J dt^J^ (^J2,,{2Vsix)J2u+i{'2Vsx) - Vij2u+l{'2\^x)J2u{2^/s^ 



X X 



yKy)'T:Preal{x) , 



(3.17) 



where in the last step we substituted s = and t = v"^ . Compared to the Hermitian limit of the 
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chiral complex ensemble the density is not proportional to 5{y) but to —ySiy)'. This reflects the 
repulsion from the real axis. Still it gives unity when integrated over the imaginary part y, after an 
integration by parts. 




Figure 3: The microscopic spectral density of complex eigenvalues at weak non-Hermiticity parameter 
a = 1.2 (left) and a = 2.5 (right), both for = 0. 

To illustrate our results further we also give the complex density for larger values of the non- 
Hermiticity parameter a. As it is seen in Fig. |21the oscillations are rapidly damped, with the density 
starting a to develop a plateau as at strong non-Hermiticity (see Fig. [TJ. For the latter as well as 
for an analytic check of this relation by letting a ^ cxd we refer to section ^ We also remind the 
reader that the overall height of the density varies with a. This is because of the normalisation to a 
delta- function in the limit a — > eq. (|3.15|) . 

3.2 Massive correlation functions at weak non-Hermiticity 

After discussing in detail the quenched spectral density we now turn to the presence of massive flavours. 
We present the spectral density for one pair of degenerate massive flavours ANj = 4 explicitly, as well as 
some characteristic polynomials or massive partition functions as examples. Higher order correlation 
functions and more flavours can then easily be obtained from the given building blocks eq. (|3.11|1 . 
following the same procedure. 

At finite- the simplest spectral density with non-zero flavour content has ANj = 4 masses. It is 
given by 

using eq. (|2.2UI) . Inserting the explicit expressions eqs. (|2.11|) and (|2.1U|) we obtain 

|2 



(z) = [z - z )w{z,z ■,t)[kn+i{z,z ) — . (3.19) 



r{4) 

\ KN^i[im,—im* 

While the first term gives exactly the quenched spectral density the second term is a correction term, 
just as for zero chemical potential When taking the large- A^ limit we have to rescale the masses 
in the same way as the complex eigenvalues in eq. (|3.2|) . 

V2 Nmf = 7]f , f = l,...,Nf . (3.20) 
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As we have mentioned at the end of subsection. \2.2\ eq. (|3.19|) includes a hmit when setting the masses 
to be real, ruf = m*j. This is because of the level repulsions both from the real and imaginary axis: 
the denominator as well as the numerator of the second term in eq. (|3.19)) vanishes. The same feature 
also holds for more general correlation functions eq. ()2.2U() . We illustrate this by Taylor expanding 
the limiting pre- kernel K^eakiiv^ ~^^*) the denominator in the imaginary part of the mass e — > 0, 
setting rj = Tjx + ie, 

liniKu,eafc(ir/^ - e, -zr/^ - e) = —N*"'^'^-, — 

e^O a [Tjx) 

+ 0(e2) . (3.21) 

The imaginary argument rotates J- to /-Bessel functions, where we could have obtained the same 
result from eq. ()3.16|) . setting x = ir] and y = ie there. Making the same expansion for the other 
terms, Kweakiz^^iri^*^) and its conjugates in the numerator, we arrive at the following expression for 
the microscopic limit of the density eq. ()3.19|) : 

/'i'L(0 = pL°L(0 - Apil,(0 , (3.22) 
with the first term from eq. ()3.13|) and the correction term to the quenched density being given by 



(3.23) 

j^ds dt^- e-2-(i+*)"'(J2^(2Vit Ohu+ii'i.fs ij) - Vtl2u+i{2V^ v)J2v{2^s C))) 
^'ds^';^ e-2^(^+*)"'(J2.(2^/itr)/2.(2V^r/) -l2.(2^/itr?)J2.(2Vir)^ - {i ^ C 
j^ds dt^j e-2^(i+*)°' (Vil2u+i{2^t r])l2ui2^s r/) - hui^^t ri)l2u+i{2^s r?) 



Here we have again called ijx = r] which is now real. Before discussing this main result of the subsection 
we also briefly give 2 examples for massive partition functions. The simplest case with 2Nf = 2 flavours 
follows from the even skew-orthogonal polynomial at imaginary mass z = im, comparing eq. (|2.14j) 
and (f^ : 

4''^^='^(m)/4°) = \mf''q2N{im) . (3.24) 

We can set the mass to be real here without taking limits, as q2N{i'nn) itself does not display any 
level-repulsion from the axis. We obtain the remarkably simple expression in the weak limit, 

z'^Liv) - j['^e--^/2.(2V^,), (3.25) 

where we have rotated the result eq. ()3.8|) from the previous subsection to an imaginary argument, 
skipping all constants (including N). As a second example we give the partition function for 4:Nf = 4 
using our formula eq. 1)2. 17() . 

z(^VrnVz(°) - i^N{im,-im*) 
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Upon using the same expansion as for the density eq. (|3.21j) we immediately obtain 

~ ^ fdt y'd^yi e-2^(i+*)"' [Vtl2u+i{2V7tv)l2u{2V^r,) - /2.(2V^r/)l2.+i(2^7?)) , 

(3.27) 

where we have again suppressed all constants. Due to the presence of conjugate pairs of complex 
eigenvalues the dependence on a cannot be factored out, similar to the same situation in the unitary 
ensembles |27| \12 \ ITH], Such a simplification only occurs in the absence of conjugate flavours, in the 
QCD matrix model with a sign problem [311 1151 ITB] . where the partition function can be normalised 
to be /i- independent. As a check both expressions eqs. (|3.25|) and l|3.27j) reduce to the know results 
jSmn] in the limit a ^ 0. 





Figure 5: The massive microscopic spectral density of real eigenvalues for masses r] = 8.74 (left) and 
r] = 4.26 (right), both at u = 0. 

Let us now discuss the massive density eq. (|3.22|) . In Fig. 0]we plot the microscopic, massive 
spectral density for two different values of masses rj = 8.74 and r] = 4.26, both at weak non-Hermiticity 
parameter a = 0.4 and = 0, and compare it the corresponding expression for real eigenvalues ^5] . 
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We have deliberately chosen the same values for the mass parameter as there, where numerical data 
for lattice simulations with dynamical fermions at = '33' were compared. Similar to having v = 2 
exact zero eigenvalues in Fig. ^the eigenvalues are pushed further away from the origin through the 
masses. Below we compare the corresponding expression for real eigenvalues at the same values of 
the mass parameters i25_. As for the quenched case the correlation functions keep their oscillatory 
structure, with the maxima located at the same positions as for real eigenvalues, at the given value of 
Q. Furthermore we observe that for larger masses the density moves toward the quenched density in 




Figure 6: The massive microscopic spectral density at weak non-Hermiticity parameter at a = 0.4 
with mass rj = 1 (left) and at a = 2.5 with r] = 4.26 (right), both at i/ = 0. 

Fig. as the infinite mass limit quenches the fermion determinant. On the other hand smaller masses 
moving towards the origin push the eigenvalues further out. In the zero-mass limit the density should 
approach the quenched density at topological charge = 2 approximately, as can be seen in Fig. [S] 
left. Although for non-zero chemical potential there is no more flavour-topology duality [3^ because 
of the explicit z^-dependence in the ii'-Bessel functions in the weight, this duality holds approximately, 
when the iC-Bessel function reaches its asymptotic, z^-independent form. For those values of \^'^\/{2a'^) 
in Fig. IHlleft, where the density is non-vanishing, this is indeed the case. 

As a final example we display the massive density at a larger value of a = 2.5 in Fig. 1^1 right. As 
seen already for the quenched density the oscillations are washed out (compared to Fig. fright). A 
further feature can be seen here. For the chosen value of a and r] the ^-value where the mass factor 
IC^ + vanishes, ib9m(^) = rj, now lies on the imaginary axis inside the shown picture. It can be 
seen very clearly that the eigenvalues are pushed away from the very location of the mass, in addition 
to the level repulsion from the y-axis. 

In the previous subsection we performed an analytical check by taking the Hermitian limit a ^ 
of the spectral density and comparing it to the known expression for the density of real eigenvalues. 
We have performed the same check here where we refer to [23 for the corresponding expression for 
the massive spectral density of real eigenvalues. Without going through the details we just comment 
of the structure to be obtained. 

Since our example eq. 1)3. 18(1 involves the two-point function it contains all four elements of the 
matrix kernel eq. (|2.8j) . the pre-kernel at different (complex conjugated arguments). In contrast 
to that for real eigenvalues the two-point or higher correlation functions are given by a quaternion 
determinant of a two by two matrix kernel with different matrix elements: a pre-kernel which is 
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basically given by eq. (|3.14|1 . its derivative and its integral (see jHO])- How can we obtain these three 
different matrix elements from one pre-kernel? Taking the limit a — > invokes a further derivative 
with respect to the imaginary part y = 9m(^), in addition to the expansion in eq. (|.S.23j) due to 
real masses. As a result we will get the real limit of the pre-kernel eq. ()3.11|) . its first and second 
derivative, correctly reproducing all terms in |^ for the real massive density. 



4 The strong non-Hermiticity limit 

The regime of strong non-Hermiticity in defined by taking the limit N ^ oo while keeping n G (0, 1] 
fixed. The rescaling of the eigenvalues is modified compared to eq. ()3.2() . keeping 

\/iV (3f?e z + i^m z) = Vn z = ^ (4.1) 

finite^. Hence the definition of the large- correlation functions is also altered, reading 



N^ooN VViV VN 



PstrongiCl,---,Ck) = Jim ^ i?Af ( . . . , ^= ) . (4.2) 



The difference in scaling can be easily understood as follows. At weak non-Hermiticity the eigenvalues 
remain close to the real axis, with their imaginary part being small compared to the real part. In 
order to obtain a macroscopic density at the origin of the order Pmacro{0) ~ A^ their average distance 
has to be 0{1/N). In contrast to that at strong non-Hermiticity the eigenvalues spread out in the 
complex plane, having a comparable real and imaginary part. To form again a constant macroscopic 
density Pmacroi^) ~ A" the average distance between nearest neighbours has to be 0(1/A^^). Here we 
measure distance with respect to the two-dimensional metric. 

A clear separation of scales between the two different limits can be defined as follows. At weak 
non-Hermiticity the microscopic density describes the universal fall off to zero of the eigenvalues 
away from the imaginary axis. In direction of the real axis along the maxima we obtain as usual 
limsR 

ez^oo Pweak 

(z) = Pmacroi^) ^ ■ In Contrast to that at strong non-Hermiticity the microscopic density 
attains the constant value of the macroscopic one in all directions, \vca.\z\-*oo Pstrong{z) = PmacroiO), 
as can be seen below in Figs. [3 and |H1 This is of course except along the real and imaginary axis, 
due to the repulsion from the Jacobian in eq. 1)2. 4() . The fall-off of the macroscopic density to zero 
in any direction (which is not seen in the microscopic regime here) can then be expected to be given 
by another microscopic function at the edge, in analogy to the Airy- kernel for real ensembles at the 
edge of the spectrum. When numerically solving the matrix model at finite but large-A^ care has to 
be taken to ensure a clear separation between micro- and macroscopic density when comparing to the 
analytic formulas. 

The computation of correlation functions at strong non-Hermiticity is more involved, as is already 
the case in the non-chiral ensemble |19j. In fact we will follow the same strategy as there, by first 
deriving a differential equation for the large- A^ pre-kernel at maximal non-Hermiticity = 1. A 
solution of this equation can be obtained either directly or by taking the limit a ^ oo of the pre- 
kernel at weak non-Hermiticity, keeping ^weak/oi fixed. Since the two methods match it is suggestive 
to recover the full range of strong non-Hermiticity p G (0, 1] by inserting a /i-dependent parameter 
through — > ^5](/i). However, we can also argue that such a rescaling can always be absorbed by 

^We comment on a ^-dependent rescaling later. 

^Because of the repulsion of eigenvalues from the real axis for /? = 4 we of course have to take this limit away from 
the axis, along the maxima of the density, see Fig. 
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redefining the scaling limit eq. (|4.1j) . Details of the derivation of the differential equation and its 
solution are given in the appendix IbI 

We begin by taking the limit of maximal non-Hermiticity ^ 1 on the skew orthogonal polyno- 
mials eqs. H'2.23() and (|2.24|) . In this limit only the highest powers survive and they become monic or 
a sum over monomials, respectively. We thus obtain the following expression for the pre-kernel eq. 
(j2.9j) in the microscopic large- limit eq. 1)4. 2|) : 

'c.irongl^.c, - ^22.+3Z.Z>r(2A: + 2z. + 2)(2A; + l)! j\T{j + + 1) ' ^ ' 

It can be shown that the function ngtrongiCi C*) obeys the following set of inhomogeneous second order 
differential equations 



1 

2 

1 



(, r*\\ N^-'^'hu (go .... 

Kstron9(4,C )U=1 = ^ (ICP^' 



where we refer to the appendix ^ for details. Due to the antisymmetry of the pre-kernel the right 
hand side changes sign. The solution of this equation can be constructed as explained in the appendix 
IHl The result is given by 

(ec)\ . - ^^_e^«^+c*^) 

poo rq 

X dq dpe-'i'-P' [J2u{2qC)J2ui2pO " J2u{2qOJ2u{2pC)] ■ (4.5) 
Jo Jo 

In order to obtain all correlation functions from eq. 1)4. 5|) we next give the weight function eq. ()2.21|) 
in the limit of strong non-Hermiticity eq. 1)4. 1|) : 



,2 



^(e^+r^) 



, /iG (0,1] . (4.6) 



Here we display the weight for any value of the strong non-Hermiticity /x. At = 1 the expression 
simplifies, the exponential factor cancels with the kernel and the argument inside the ET-Bessel function 
reduces to unity times |gp. Inserting this into eq. (|2.in|) we obtain the following result for the 
microscopic spectral density at maximal non-Hermiticity 



Pstron,m,=i = ^(r^-c^)iepi^2.(ici')e^^^'^^"^ (4.7) 

poo rq 

X dq dpe-'i'-P'[J2u{2qnJ2u{2pC) - J2u {2qO J2u {2pe)\ ■ 

Jo Jo 



As an analytic check we have performed the limit a ^ oo of the microscopic spectral density at weak 
non-Hermiticity eq. (|3.13p . where we refer to the end of appendix^for details. We obtain exactly the 
same functional form of the spectral density, eq. (|B.20() (up to a different constant normalisation), in 
terms of the variable = z / [i^/2. 

Together with eq. (|4.6|) this leads us to conjecture that the //-dependence at strong non-Hermiticity 
with [i < \ can be restored as follows, replacing 

e ^m^e (4.8) 
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Figure 7: The microscopic spectral density at maximal non-Hermiticity /U = 1 with v = Q (left) and 
V = 2 (right) . 



everywhere in eq. (|4.7() . We have shown it to be true for p = 1 and for small p <^ 1. However, if this 
is true we could then redefine our strong non-Hermiticity limit eq. (|4.H) to completely reabsorb this 
//-dependence, defining 



p' 



pV2 



^/N z = ^ 



(4.9) 



ik ,,2k 



Pstrongi^lj • • • ^ ^k) — 



2V- 

hm KT-, — r 

N~.oo (1 + fi^fN^ 
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iiP^/2 



ikpV2 



As an illustration of our results we show the microscopic spectral density eq. 1)4. 7|) and its de- 
pendence on V at maximal non-Hermiticity // = 1 in Fig. Q One can clearly see that for v = 2 the 
additional zero eigenvalues push the density away from the origin. For the numerically evaluation 
of the integral in eq. 1)4. 7(1 necessary to plot it an equivalent representation of the integral is very 
convenient. Using the identity derived in the appendix eq. (|B.23|) we have 



PstrongiO — ('^ 



* 2 



dr 



_/2.(riensinh(iyr 



*2- 



(4.10) 



where the exponential pre-factor has cancelled out. This form of the density allows to dicuss its 
rotation and reflection symmetries most clearly. The fact that it is an even function of (^^ — C*"^) ^^'^ 
otherwise only depends on the modulus leads to the following properties. First, one can obviously 
rotate the density by it/2, vr, or 37r/2, multiplying ^ by i , (—1), or —i without changing it. Second, 
the density is invariant under reflections along the x-, y-, (x = y)- and (x = — ?/)-axis where ^ = x + iy. 
For example for the (x = y)-axis in the first quadrant we write the pair of ^-values related by reflection 
as 

^± = re'^^^i^ ^ C± = -i^T- (4-11) 

The modulus remains unchanged, |^+| = |^_|, as well as the difference (^4. — ^!!J_^) = —{C-'^~^~)- This 
symmetry can be observed in Fig. |H1 The remaining reflection symmetries easily follow in a similar 
way. 
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Figure 8: Part of Fig. [3 (left) showing the first quadrant only. The symmetry with respect to the 
diagonal SRe.^ = Qm,^ is clearly visible . 



4.1 Massive correlation functions at strong non-Hermiticity 

In this subsection we give the simplest nontrivial result including one pair of degenerate massive 
flavours, 4:Nf = 4, where we will follow closely the corresponding subsection \'A.2\ at weak non- 
Hermiticity. Higher order correlations as well as more flavours follow along the same lines. We 
start by defining how to rescale the masses 



Nruf = r]f , f = l,...,Nf, (4.12) 

in analogy to the eigenvalues in eq. 1)4. 1|) . In order to derive the massive spectral density we could 
simply take the limit a ^ oo of the expressions in eqs. (|3.13|) and (|3.23j) at weak non-Hermiticity, 
together with the results from appendix IbI 

In order to give more compact expressions we derive results using the alternative form of the strong 
density eq. (|4.1()|) and the corresponding kernel. The massive spectral density follows again from eq. 
()3.19|) at finite- A^. The kernel at strong non-Hermiticity we use is obtained from eq. ()4.5|) by inserting 

(e,r) = / ^^^^/2.(.e.a)sinh(-v^(el-a^)j • (4.13) 



eq. (|B2 



/^strong 

Expanding again in the imaginary part of the mass e — > 0, setting rj = r]x + ie we easily obtain for the 
denominator of eq. (|3.19|1 



liuiKstrongiiVx - e,-i'nx - e) = ^ ^\^^ {-2i) € r]^ J drhvirril) + 0{e^) . (4.14) 



The same expansion can be done for K-strong 

(^,±ry(*)) and its conjugates, leading to the following 
expression for the massive spectral density at strong non-Hermiticity: 

PstlongiO — pitrongiO ~ ^P^trongiO ) (4-15) 
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with the quenched density given by eq. (|4.1fl|) . The correction term reads 
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strong 
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4tt 



dr 



fJ2u{r^'n) sinh[^\/l - (^^ -\- ■q'^')] 



/ ^=^J2u+i{rCv) sinh[-Vl-r2(r2 + r?^)] - drT^J^.^rCri) cosh[-^l - r^{C^ + 
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7] I dr huirrj^) 




(4.16) 



We have again cahed r] = rj^ which is now real. Before discussing the massive spectral density let us 
also give some examples for partition functions. Taking the strong non-Hermiticity limit of eq. 1)3. 26() 
at 4:Nf = 4 we simply have to insert the expansion eq. (|4.14|) to obtain 



strong \ IJ 



~ / drl2u{rrj^) , 




(4.17) 



up to suitable normalisation constants. Since we did not compute the asymptotic even skew-orthogonal 
polynomials at generic strong non-Hermiticity < /i < 1 we cannot directly use eq. (|3.24|) for 2Nf = 2. 
We can circumvent this problem by taking the infinite a limit on eq. 1)3.25(1 . leading to the following 
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strong 
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Let us come back to the spectral density eq. ()4.16() . Below we show a few examples for the influence 
of the mass parameter at maximally strong non-Hermiticity. Compared to the quenched density in 



P strong ('^) 




P strong (0 



4 0.1 



. 05 



9m ^ 




9m ^ 



Figure 9: The massive spectral density at at mass r/ = 2.5 with z/ = (left) and u = 2 (right) . 

Fig. [7|left there is an additional level repulsion from the mass, when ib3'm(^) = 77. Non-zero topology 
leads to an additional level repulsion from the origin as can be seen in Fig. right, and both the 
effects of mass and v can be clearly separated (see also Fig. [7| right). 

For smaller values of the mass parameter the density approaches approximatively the quenched 
density at v = 2, as the mass acts as additional zero eigenvalues (comparing Fig. ^Jleft and Fig. [T] 
right). On the other hand going to larger masses we get back the quenched density as can be seen in 
Fig. [TUl right. 
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Figure 10: The massive spectral density at at mass rj = 0.5 (left) and r] = 8 (right), both with u = . 

Of course there is still level repulsion at the large mass value, in this case r/ = 8, as displayed 
in Fig. 111! However, this does no longer affect the density at the origin. Furthermore, there are 
no strong oscillations observed here at the location of the mass, as seen in unquenched QCD jl6j . 
Here, the situation is very similar to phase quenched QCD also reported in ^] on the matrix model 
in the (3 = 2 symmetry class. These findings underline once more the importance of the presence or 
absence of the sign problem on the level of the microscopic complex eigenvalue correlations of the Dirac 
operator. We note in passing that because of the analytic relationship between the correlations at 




Figure 11: Part of Fig. (right) showing the first quadrant for a larger argument 9m(^). The level 
repulsion at the mass value ^ = 8z is clearly visible. 

strong non-Hermiticity and weak non-Hermiticity at large a similar figures could have obtained in this 
section by plotting the corresponding pictures at large enough a. We have checked this numerically 
in several cases for both the quenched and massive density, with a = 5 (see also Fig. ISJ. However, 
such a density at weak non-Hermiticity will always vanish for large enough ib3'm(^), at any given a. 
In contrast to that the density at strong non-Hermiticity remains constant. 
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5 Relation to the Dirac operator spectrum with chemical potential 



The purpose of this section is twofold. First, we give a more detailed discussion of the matrix repre- 
sentation of our complex eigenvalue model eq. (|2.4j) in terms of two non-Hermitian quaternion real 
matrices, including mass terms. This justifies to speak of a matrix model ensemble. It also gives an 
explicit realization of the joint probability distribution in eq. (|2.4() by computing the Jacobian that 
leads to complex eigenvalues. In particular it yields the special weight function eq. (|2.21j) containing 
a i^-Bessel function. The second purpose is to relate the model to the Dirac operator spectrum of 
QCD-like gauge theories in the presence of a chemical potential. To do this matching we need to 
compare the global symmetries of such Dirac operators to a random matrix representation. 

The Dirac operator as it appears in four dimensional gauge theories can be represented by an anti- 

Hermitian block matrix, Ipgauge = ^ ^ ; with non vanishing entries only on the off diagonal 

due to chiral symmetry. For different gauge groups and representations the off-diagonal blocks W have 
either a orthogonal, unitary, or symplectic symmetry, as was pointed out in jH]. Adding a chemical 
potential /U for the quarks amounts to shifting the Dirac operator by Ipgauge — > Ipgauge + /^7o where 
7o is a Dirac-7 matrix. In a given basis it is represented as a block matrix as well, with unit elements 
on the off diagonals only Being Hermitian this addition destroys the anti-Hermiticity property, 
where the sum Ij) gauge + A*7o now has complex eigenvalues. The generalisation of the classification [HI 
to /i 7^ we use here was made in jlO| . 

If we replace the Dirac operator by a random matrix we have to satisfy two criteria. First, Ipgauge 
may have exact zero eigenvalues corresponding to different sectors of topological charge through an 
index theorem. For a fixed number v of zero modes, which we shall always consider here, we thus have 
to replace the blocks W by rectangular matrices, of size {N + v)xN. Second and more importantly, 
these rectangular matrices have to be in the same symmetry class as the corresponding gauge theory, 
consisting of real, complex or quaternion real elements [01 E3]- The latter case which we consider here, 
also labelled by the Dyson index /? = 4, describes SU{Nc) gauge theories in the adjoint representation. 
When representing the gauge theory on a lattice using staggered fermions the corresponding symmetry 
is an SU{Nc = 2) gauge theory in the fundamental representation. 

The addition of a chemical potential in the matrix model can be done in different ways. It was 
first suggested to add it as a constant times 70, the off diagonal unity matrix, assuming this 

term to be diagonal in the random matrix space. While this model successfully explained the failure 
of the quenched approximation an analytical treatment of all complex eigenvalue correlations 
from orthogonal polynomials similar to the n = case was not achieved to date, apart from the 
quenched spectral density obtained using the replica method jl2j . Therefore a different matrix model 
representation of the chemical potential term was suggested very recently in ^21 for the QCD symmetry 
class (5 = 2. It assumes that the chemical potential term is non-diagonal in matrix space and is 
represented by a second, uncorrelated matrix which has the same symmetries as the one modelling 
the /i = part. Making the model more complicated at first sight it has the remarkable feature of 
leading to a complex eigenvalue representation, which then allows for a treatment using (bi-)orthogonal 
polynomials |15|ll6j. A different model given only in terms of complex eigenvalues was first suggested in 

and tested in quenched QCD lattice simulations jl3j . Although it is lacking the matrix symmetries 
of Ipgauge it asymptotically agrees with the results of ^lEI (see the weight eq. (|2.22|) vs. eq. (|2.21j) ). 

In the following we will adopt the strategy of and apply it to the complex /3 = 4 symmetry 
class. Our corresponding matrix model is thus given in terms of two rectangular (A^ + z^) x N matrices. 
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$ and ^, with quaternion real elements without further symmetry properties: 

^'''hi^f}) - / ^^rf^exp [-iVTr(cI>tci> + vl;t^)] JJ det ( ^^(^f^^^ '^^^f j , (5.1) 

where 1 is the quaternion unity element. Compared to eq. (|2.ip we have now included the quark 
masses. As a first step we define the linear combinations 

C = +[1^^ 

D = i<^^ + /i^''^ , (5.2) 
with a trivial constant Jacobian. In terms of these new matrices the partition function reads 

.Aril±^Tr(CtC + DD^) - N ^^~ f^ iiDC + C^D^) 



X 

/= 

In this form we have only convergent integrals for /j, £ [0, 1] at finite- A^, as can be seen from eq. 1)5. 8p 
below. We will not discuss the possibility of phase transitions here, where for lattice simulations of 
this situation we refer to |35j . 

The Dirac matrix inside the determinant has i' zero eigenvalues and complex eigenvalues which 
come in pairs of opposite sign as well as complex conjugate pairs. The former is due to chiral symmetry 
and the latter is due to the symplectic structure as the eigenvalues of a quaternion real matrix come 
in complex conjugate pairs. It is this symplectic structure that is responsible for the fact that the 
determinant will remain positive definite, without having a sign problem as unquenched QCD. It 
makes this symmetry class an ideal testing ground for lattice simulations with dynamical fermions in 
the presence of a chemical potential. 

In the next step we parametrise the matrices C and D as follows, 

C = U{X + R)V 

D = V\Y + S)U^, (5.4) 

which is equivalent to making an independent Schur decomposition oi CD and DC. Here U and V 
are symplectic matrices. The matrices R and S are upper triangular with real quaternion elements, 
and we have split off their diagonal parts X and Y which are also quaternion real, without containing 
the second and third quaternion unit (see e.g. J7j for explicit representations). They contain the 
complex eigenvalues of C and D which come in N pairs (xfc,x^) for the matrix X and {ykiVX) 
the matrix Y respectively. The parametrisation is not unique as under an additional diagonal unitary 
transformation the matrices remain triangular. It only changes the upper triangular matrices R and 
S which we will integrate out afterwards. This symmetry can be used to restrict the matrices to be 
V e Sp{N)/U{l)^ and U G Sp{N + u)/{Sp{u) x U{1)^). The parametrisation eq. ((Ol) as well 
as the counting of degrees of freedom is discussed most explicitly in appendix O The relation to the 
original eigenvalues of the Dirac matrix is 

XkVk = -4 , (5.5) 

which explains that the eigenvalues of CD come in complex conjugate pairs (z^,z^2), k = 1,...N. 
The final Jacobian of the transformation computed in appendix ICl reads 

i six^^ms) ) - n n \4 - 4f \4 - n \4 - 4¥ , (5.6) 
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after a suitable ordering of the independent variables. We can now integrate out the symplectic 
matrices U and V , as well as the Gaussian integrals over R and S to obtain 



Af/ iV Af 

n + p n - -'I' - T n i-^ - . (5.7) 

/=i fe>/ /i=i 



X Ixjl"^" 



Here we have included the mass dependent factor ~ \mf\'^'^ in the normalisation, to the power of 
mass degeneracy 2 in 2Nf. Although we have now achieved an eigenvalues representation in terms of 
Xk and Uk we still have twice as many degrees of freedom that what we aimed at, a single set of 
complex eigenvalues Zk- If we substitute yk by Zk'. yk = —z^/xk we can integrate out the variables Xk 
in analogy to 13", after decomposing them into polar coordinates. 



JV(1-,,') 



2 , ^*2n 



4a. 



{z^ +z 



2 ^3 ' J 



(5.8) 
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Nf 
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if=l ^ ^ / 7(0) 



J|det 




This gives us a complex eigenvalue model as in eq. (|2.4p with the specific weight function eq. (|2.21j) 
plus mass terms, as well as a matrix representation of the q2N{z) eq. 1)2. 14() for Nf = 1. We note 
that compared to the QCD symmetry class the index of the A'-Bessel function is shifted: v — > 2u. 
The massive partition functions eq. ()2.16() computed in section 12.21 are obtained by taking ANj flavours 
of twofold degenerate real masses m/ in eq. (|5.1j) . or 2Nf flavour pairs of complex conjugate masses 
(mf,m*ji) for complex masses: 

~ (n/- ( ) ( ™> ) ) . (5..) 

Let us stress that it is not this imposed degeneracy that is responsible for the absence of a phase in the 
weight. The non-degenerate partition function eq. (|5.8() has already a positive definite weight. Eq. 
(|5.9j) provides a determinant representation for the mass terms in eq. (|2.1H|) for which we have been 
able to compute all complex correlation functions. The computation of partition functions without 
2-fold degenerate masses eq. (|5.8|) remains an open problem for 2Nf > 2, where the special case 
2Nf = 2 is given by eq. ()3.24() . In the limit /x — > the partition function eq. 1)5. 9|) reduces to 

lim Z^^^f\{mf}) ~ / II dxj exp [-Nxf^ "Q |^^|4-|^2 ^ ^2 14 "Q |^2 _ ^2|4 ^ 

using eq. (|3.15j) . This is the partition function of the chGSE with real eigenvalues computed in 
|251 132j and successfully compared to lattice simulations with dynamical fermions at ^ = in [331 125| . 
It matches the corresponding chiral Lagrangian in the e-regime calculated for equal masses in |36j . 
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6 Conclusions 



In this paper we have solved the complex extension of the chiral or Laguerre Symplectic Ensemble 
to non-Hermitian matrices. The solution was shown to be expressible in terms skew orthogonal 
polynomials for general weight functions in the complex plane. We gave explicit expression for finite- 
N for all correlations functions of eigenvalues in the absence and presence of twofold degenerate mass 
terms, as well as for characteristic polynomials. Two examples of a Gaussian Ginibre type weight and 
a weight with a /C-Bessel function were given. We proved that the (skew-) orthogonal polynomials 
corresponding to the latter are Laguerre polynomials in the complex plane. 

We investigated the large- limit of complex skew orthogonal Laguerre polynomials, first treating 
weak non-Hermiticity, and derived all correlation functions in terms of the limiting kernel and weight 
function. These finding were illustrated by examples of the microscopic spectral density without and 
with mass terms as well as by partition functions, where we discussed the influence of both exact zero 
eigenvalues and the masses. The same analysis was then repeated for strong non-Hermiticity where 
we had to solve two inhomogeneous differential equations to determine the limiting kernel. The virtue 
of the weak non-Hermiticity limit is that it allows to extrapolate in terms of the weak non-Hermiticity 
parameter a between real eigenvalue correlations at a — > 0, and strong non-Hermiticity in the limit 
a — > oo. This provided several analytical checks. 

These mathematical results may have applications in various physical systems with complex op- 
erators, and we have focused on the application to the Dirac operator spectrum in QCD-like gauge 
theories. We have identified the correct global symmetries of our ensemble by finding a two-matrix 
representation that matches with the symmetries of the Dirac operator with chemical potential in 
SU (Nc) gauge theories in the adjoint representation, or SU {Nc = 2) gauge theories in the funda- 
mental representation using staggered fermions in a lattice discretisation. Our two-matrix model 
representation led to the particular form of a i^-Bessel weight function of complex eigenvalues, which 
is identical to that of the corresponding model in the QCD symmetry class. However, the solution we 
presented in terms of complex eigenvalues holds for more general classes of weight functions as well. 

It would be very interesting to verify our predictions in lattice simulations with dynamical fermions 
in the corresponding theories with chemical potential. In contrast to QCD these do not suffer from a 
sign problem and thus can be solved using standard techniques. This project is work in progress. 
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A Orthogonality of complex Laguerre polynomials 



In this appendix we prove that Laguerre polynomials in the complex plane are orthogonal with re- 
spect to the weight function eq. 



(|2.21j) . This justifies our construction of skew orthogonal Laguerre 
polynomials described in subsection 12.31 We wish to show that 



\2u+2 



Ki^{a\z\'^) exp 



{z' + z*') 



2b 



^* 2 



2b 



^kl , 



(A.l) 
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holds for a > b when integrating over the full complex plane , where we have introduced 

iV(l-/.2) 



a = — 75 and b 



2/^2 



(A.2) 



Note that throughout this appendix we use the shifted topological index 2z/ — > as it appears in the 
complex chUE [El, compared to the weight eq. (|2.21j) . In the following we can restrict ourselves to 
the case A; > Hn eq. ()A.1|) . as then k < I follows from complex conjugation. To show eq. ()A.1|) it is 
thus sufficient to prove 



(fz \z\'^''+^K^ia\z\'^)exp 



2^^ J 



2b 



, for all /c > / 



(A.3) 



The remaining case A; = / in eq. (|A.1|) is trivially non-vanishing, as we integrate over a positive 
quantity. The corresponding norm of the k-ih. polynomial is computed below as well. To set up our 
orthogonality proof a la Gram-Schmidt we first compute the following general integral. 
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(A.4) 



After introducing polar coordinates we can integrate out the angle to obtain the modified I-Bessel 
function. Changing variables = s the remaining integral can be found in terms of the hypergeometric 
function j^HI where convergence is guaranteed due to a > 6. We note that all the moments Cki are 
real. The case k < I can thus simply be obtained from interchanging k and I: 



Clk 



Ckl for I > k . 



(A.5) 



We first consider the orthogonality of the zero-th polynomial to L^. For / = the hypergeometric 
function in eq. ()A.4|) reduces to a simple rational function, F{n, (3; (3; z) = (1 — z)"", and we obtain 
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Using the fact that the Laguerre polynomials have an explicit representation, 
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we can simply write down the orthogonality relation between the k-ih. and 0-th polynomial as 
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Here we have used the definitions (|A.2|) as well as ^-^^2) = °^ 2h^ — = N"^ / fi^ . The insertion 

of the coefficient Cmo eq. ()A.6|) leads to the binomial sum (1 — 1)'^ that vanishes for all k > Q. For 
= we obtain the norm of the 0-th polynomial reading 

= ttn-'^-''{i + . (a.9) 

In the general case eq. (|A.3|) the strategy of our the proof will be similar. We shall try to reduce the 
expressions to vanishing binomial sums. However, there will be remaining terms, which we have to 
show to cancel as well. As a preliminary step we can use a transformation formula |29j to show that 
the hypergeometric functions in the coefficients c^i eq. HA.4|) terminate, 

F[l + k + v,l + k;l + k-l;-^ \ = F(-l-v,-l-l + k-l--^\[l--^\ . (A.IO) 

Because of > 0, and k > I > being integers the hypergeometric function on the right hand side 
terminates, containing / + 1 terms. Writing down the definitions we obtain for eq. (|A.3|) 

m=0 ^ ^ ^ / 

^ \k-m m\\ 2b J + Z^^ ) [k-mjm\[ 2b 

= Si + S2 , (A.ll) 

where we have split the sum into two contributions for m < I and m > / in Cmi- In the sequel we 
manipulate the Si and S2 separately by rearranging the hypergeometric functions, showing that they 
indeed cancel for k > I. We start with the simpler second sum, using eq. ()A.10|) : 
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In the next step we reduce the individuals sums sorted by powers 
binomial sums by shifting the summation indices: 



^2 = ^^f^ (1-^2). (-) I 
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In the last step we have used E1=;^+I"(-)'" = - YZ,=^{-T ^'"l^^"") . replacing incom- 

plete binomial sums. In particular we need k > I for the first term. Also we have multiplied in the 
factor 'l^_^2y ■ The generic n-th term in the sum is explicitly displayed. In this form we can compare 
to the first sum Si, again after doing some manipulations: 

m=0 \ / \ / 



(k-m)\{m^v)\m\\ 2b J (l - m)\ a'^^+'^+'' 

/ ^2x-i-m-;-. . 
X 1 ^ F —m — u, —m; 1 + I — m; 

A^l^^/^^ ^ ^^J/fe-m)!(m + z/)!m!(/-m)! (l-/i2)2m 

mim + v) ^2 J7^(j7^ _ i)(77), _|_ j/)(77i _l_ J/ _ 1) 54 m!(m + — m)! 52*' 

^ 1 + T. — ~, H 7:; — ~^ N . , tt: J + . . . + 



(l + /-m)a2 (l + Z-m)(2 + /-m)2 l\u\m\ a?"^ , 

In contrast to eq. ()A.12|) the length of each hyper geometric sum now increases from term to term. 
Let us therefore write them out most explicitly: 
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In order to see that Si and S2 cancel we now have to compare individual powers p- = j^^rj;^- For 

2{;-l) Ci I „2\2{I-1) 

the for highest power, = (i-|t2j2(i-i) > it is easy to see that this is indeed the case, with only the 
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first term in tlie last line of Ei in eq. ()A.15|) contributing. In contrast to that for the lowest power |^ 
the last term in each hypergeometric sum in eq. HA.15|) does contribute, leading to 

This clearly cancels the last term in eq. ()A.13|) . To see that all other terms cancel too we now write 
down all contributions to the generic power ),_ 2(21-2^-2 to be compared with that of E2 in eq. ()A.13|) . 

In El we thus have to pick all terms l2q~2p with 2q — 2p = 2l — 2n — 2. Because of n < / — 1 and p < q 
only those hypergeometric sums with 2q > 2/ — 2n — 2 start contributing. Eliminating p = q — l + n+l 
we obtain the following contribution 
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where we have shifted summation variables in the last step. This term clearly cancels the corresponding 
term in E2 in the last but one line of eq. (|A.13|) . This ends our proof of eq. (|A.3|) and therefore of 
the orthogonality of the Laguerre polynomials in the complex plane. 

What remains to be computed are the norms of the k-ih. Laguerre polynomial, which can be 
obtained as follows. The only step where we used k > / was when producing a zero in the first 
term of E2 in eq. (ICTll . deducing (1 - l)''"' = 0. For k = I this is no longer true and this 
term does contribute, leading to the only non-vanishing contribution from Ei + E2. Because of the 
orthogonality only L'^{z'^)z*'^'^ will contribute to ||L^|p. To obtain the norm for the polynomials in 
monic normalisation, 

we still have to normalise eq. (|A.11|) appropriately, leading to 



Nz^ 

~2 



which is our final result on orthogonal Laguerre polynomials in the complex plane. 

In order to clarify the issue of orthogonality of complex Laguerre polynomials with respect to 
different weight functions we have explicitly checked that Lq and are not orthogonal with respect 
to the exponential Gaussian weight in eq. (|2.22|) times Izp'^^^ for a general value of as was claimed 
in . The integral is given by a linear combination of Legendre functions of index v + \ and + | 
which does not vanish in general. Only in the special case of = ib^ where the Laguerre polynomials 
are related to Hermite the orthogonality holds. This does not come as a surprise now in view of the 
simplification of the i^-Bessel weight eq. (|2.2H) to eq. (|2.22j) in that case. 
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For comparison we shall also give the norms for the Gaussian weight eq. (|2.27|1 . It trivially holds 
that for all rotational invariant weight functions w{z,z*) = w{\z\) the monic powers form a set of 
orthogonal polynomials 

J cPz w{\z\)z^z*^ ~ 5m . (A.20) 

For different weight functions only the corresponding norms differ, which then enter the kernel, or the 
construction of the skew orthogonal polynomials as in section 12., 31 For the Gaussian weight eq. (|2.27|) , 



w'^^\z,z*) = Izl^^^+^expi-A^lzpJz^fc , (A.21) 
with shifted 2z^ ^ one easily obtains the following result, after changing to polar coordinates 

= /"d^z |zp-+iexp[-7V|z|2]z2fcz*2fc ^ ^ r(2fc + t/ + |) ^^^^^^ 

B Differential equation for the kernel at strong non-Hermit icity 
B.l Construction of the differential equation 

In order to derive the differential equation we first rewrite the coefficients inside the pre-kernel eq. 
()4.3|) in a more symmetric way. Using the following identity for the F-function called doubling formula 

T{z) (2^)1 



T{2z) 22--^F(z + i) ' 

we obtain 



(B.l) 



^stron,(e,r)l^=i = 24.4-4^ EE 23(fc+.)(2A: + 1) ! ! (2j) !!F(A; + \ + v + l)F(j + + 1) ' ^^'^^ 
It is convenient to change variables, 

V = (B.3) 

and to define the following functional which is antisymmetric in the variables n, v, 



oo k 

aiu,v) = EE 

k=0 j=0 

This leads to 



2fcF(A: + 1/ + |)(2fc + l)!!2ir(j + i/ + l)(2i)!! 



(B.4) 



jY2i^+2 

Kstrong{CX*)\^i=l = 2i'^+3^ (t{u,v) . (B.5) 

Applying the operator v~'^'^ d^v'^'^^^ to a{u,v) relates to va{u,v). In detail we have 

°° '^-i ,,2fc+l „,2i+l 

\v-^''dy+^dJa{u,v) = yy-, (B.6) 

^ -J V , ; A.A^2feF(A; + z.+ |)(2fc + l)!!2Jr(j + i. + l)(2j)!! ^ ^ 



oo k 



EE 



y2k y2j 



2*^-iF(A; + iy + l){2k- 1)!! 2ir(j + u + l)(2j)!! 



fe=0 j=Q --^ V" i - I 2 
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where we have shifted the sum over j in the first term and denote (—1)!! = 1- This compares to 
multiphcation by v reading as follows 



V a{u, v) 



u 



2k+l 



2^T{k + v + ^)(2k + 1)!! 2ir(j + u+ l)(2j)! 



oo k 

Y.T. 

k=Oj=0 
oo fc—l 

~ § ^ 2^r(A; + + i)(2A; - 1)!! 2JT(j + u + l)(2j)!! 



,2i 



(B.7) 



where we have shifted the sum over k in the second term. Subtracting this from eq. HB.6|) we obtain 

\2k+l 



[uv 



2^T{k + U + |)r(A; + + 1)(2A: + 1)! 



2(to)° 



r(z. + i)r(z. + i)"^22'=-ir 



2fc 



+ z^ + i)r(A: + i/ + l)(2A:)! 



22i/+l 



^ r(fc + 2i/ + 



l2v{2y/m)) 



vr (to)'' 



fB. 



In the first step we have again used the relation (|B.1|) and the last step follows upon the series 
representation of the /-Bessel function. Going back to the kernel with eq. HB.5|) in the changed 
variables eq. ()B.3|) the differential equation (|4.4|1 follows, 



l^strongiCjC )\fi=l 



N^'^+'^hu (CO 
4vr i^C*)^" 



(B.9) 



c-2 

The corresponding equation in the variable u = ^ is obtained in the same way, or by simply using 
the antisymmetry of the function a{u, v) 



'^strong 



+ 



2u+2 



4vr (CC 



*\2v 



(B.IO) 



B.2 Solution of the differential equation 

Next we seek for a solution of these two inhomogeneous differential equations. The homogeneous 
equation can be easily solved, as it holds 



, 



(B.ll) 



in both variables u and v. The second independent solution in terms of i^-Bessel functions is excluded 
because of the regularity of the pre-kernel at the origin. We therefore conclude that the following 
product 







Iy{u)Iy{v) 

{uvY 



(B.12) 



33 



solves both homogeneous equations. However, this solution is symmetric in both arguments, in contrast 
to the antisymmetric pre-kernel. Each factor of the homogeneous solution has the following integral 
representation that will be useful below. In terms of the original variables it is given by eq. (|B.3|1 



I dqe-"'^J2u{2qO 



(B.13) 



and similarly for d^* . 

Instead of making an Ansatz that solves the full set of inhomogeneous differential equations ()B.9|) 
and HB.10|) we construct a solution by rewriting the right hand side as an integral and manipulating 
it until we obtain the correct differential operator acting on it. Omitting the constant pre-factor we 
start from the following integral identity. 



hu (CO 



4e^(«'+f*') 



/ dqe~^'^ qJ2u{2qOJ2u{2qC) 
) Jo 



_4e^(C^+C-)l 
2 



(B.14) 



dge-'?' [(1 - 2q^ + qdg)J2,{2qO] / dp e~P' J2,{2pC) 



*\2u 



(4. + l),^_^2 



2 ^ V" ; ^^ d. 



^2 



-4e^(«'+f*') 



dqe-1 J2u{2qi) / dp^-^ J2u{2pC) 



(ec) 



*\2u 



dqe-'i J2ui2qC) / dpe"" J2u{2pC 



Here we first integrated by parts e J2i/{2q(*), with a vanishing surface term. Next we replace 
qdqJ2u{2q^) = ^d^J2u{2q^), and then use that J2u{2q^) satisfies a Bessel differential equation: 
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J2ui2qO 



-2q^J2u{2qO ■ 



(B.15) 



We can now take the differential operator out of the integral and interchange it with the ^-dependent 
pre-factor, giving us precisely the desired differential operator. This provides a solution of the sec- 
ond differential equation (|B.10|) . after multiplying with the appropriate constant. A solution of the 
corresponding equation ()B.9|) for (* follows along the same lines 



(ec) 



5, 



C 



* 2 



e-" J2.(2gO / dpe-P J2.(2pe) 



^ (ec*)2- 

(B.ie) 

However, none of two solutions is antisymmetric in the variables ^ and C* to be proportional to the 
pre-kernel, nor do they obviously solve both differential equations with the different signs simulta- 
neously. We will therefore add and subtract a multiple of the homogeneous equation eq. HB.12|) in 



its integral representation eq. (|B.13|) to eq. ()B.14|) and from ()B.16|) . respectively, in order to obtain 
an antisymmetric solution of both differential equations. The requirement of antisymmetry of the 
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pre-kernel fixes the solution uniquely. It can be seen that the solution is given as follows 



/ dq / dpe-'^"'P" [J2u{2qC)J2u{2pO " J2ui2qOJ2u{2pC)] (B.17) 
JO Jo 



poo pq 

/ dq / dpe-'^'-P' J2u{2qOJ2u{2pC) 
Jo Jo 



poo pq 

/ dq / dpe-'^'-P' J2u{2qQ*)J2u{2pO 
Jo Jo 



l^y^ dqe-'^'j2.{2qC)\r-^J^ dp^-^\hA2vC 



The first line is antisymmetric in its arguments ^ and In the subsequent lines we have made an 

2 2 

integration by parts of either e"*? J2u{2qC*) or — e""^ J2v{2q^) in order to show, that this expression 
solves both differential equations, ()B.14|) and ()B.16|) plus or minus a solution of the homogeneous 
equation. The full solution for the pre-kernel including all factors thus reads 

]\j2u+2 1 



2tt (^C*)^^ 

poo rq 

/ dq dp e-'i^-P' [J2,{2qQ*)J2,{2pi) - J2u{2qi)J2u{2pQ*)] . 
Jo Jo 



As a check and alternative way to derive this result we take the limit a ^ oo of the pre-kernel 
at weak non-Hermiticity eq. (|3.11|) . Substituting first t hy p^ = 2sta^ and then s by = 2sa^ we 
obtain from eq. ()3.11|1 

hm K^eafc(e,r) = lim ^^^''^^'tS^; (B.19) 

a— >oo a^oo Or j 



tV2 

dq I dp e-«^-P^ 



(2qQ*\ , ( 2pi\ ( 2qi\ ( 2pC 

j j " I WlJ \^2 

If we then replace the first integral dq — > dq, keeping and fixed, we reproduce 

the strong kernel eq. ()B.18|) up to constants and the exponential pre-factor. Therefore we cannot 
directly identify the kernels, but have to take into account the weight function as well when mapping 
correlation functions from the weak to the strong limit. From eq. (|3.12|) we thus obtain 

hm 2ivp^,,fc(0 = T^(s'-ei)ie5i'i^2.(ie5i')e^(«i+«") (B.20) 



/"oo pq 

/ dq dp e-^'-f' [J2. {2qes) J2. (2^5) - (2^?!) ^2. {2q^s)] , 
Jo Jo 



exactly reproducing the strong density eq. (|4.7|) obtained from eq. (|B.18|) . Here we have introduced 
the scaling variable 



with ^ being the scaling variable in the weak limit eq. ()3.2() . ^5 scales as the scaling variable in the 
strong limit eq. ()4.ip . The explicit factor 2N in front of the density is added since the weak and strong 
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density are rescaled with different powers of N, see eqs. (|3.8|) and (|4.2j) . The exphcit /i-dependent 
factor in front can be absorbed by defining a /^-dependent strong scahng Umit by ^5, and rescaUng 
the density accordingly. It can be easily seen that all higher eigenvalue correlation functions at weak 
and strong non-Hermiticity can be matched in the same way when taking a ^ 00. 

We end this appendix by noting that another integral representation can be obtained from the 
a ^ 00 limit of the weak pre- kernel eq. First we substitute s = and t = v^. Next we 

interchange the integrals with respect to u and v and substitute uhy = v?a'^{\ + v^) to arrive at 



lim 2a^ I ds 



1 (k_ 



-2s(l+i)a2 



lim 4 



dv- 



1 + 



dv- 



1 + 



i2u 



dy y 



J: 



2yvi 



2u 



J: 



2yC 



2u 



ayj2{l + v'^) 



a 



2(1 



(B.22) 



In the last step we have replaced j^^'^^^'^'" ) _^ which case the integral can be evalu- 

ated. Keeping the same ratios -^^^ and fixed, we obtain the following useful integral identity by 
comparing eqs. ()B.19|) and (|B.22|) : 

poo pq 

4 / dq dp e-^'-f' [J2u {2qCs) J2. (2^5) - ^2. {2qCs) ^2. (2K5)] 
Jo Jo 



dv- 



I2u 



2 e-^«l+CS^) 



dr 



2vCsC:s 

1 



(l+l)^) 



(l+li^) 



/2.(re5C5)smh Qyr37^(e| - cf)^ . (B.23) 



In the last step we have substituted r = 2v/{l + v'^). It is this last form which is very convenient for 
a numerical evaluation in order to plot the spectral density at strong non-Hermiticity. 



C Jacobian of the two-matrix model 

In this appendix we compute the Jacobian for the parametrisation eq. (|5.4|) 

C = UTV , T = X + R 

D = V^QU^ , Q = Y + S , (C.l) 

resulting from Schur decompositions of CD and DC. To this aim we proceed in two steps. First, 
we provide a labelling of the independent degrees of freedom to obtain a block triangular form of the 
Jacobi matrix in terms of quaternion matrix elements. Here we follow some analogy to T7'| appendix 
A. 35 where a similar decomposition was given for a single quadratic complex non-Hermitian matrix. 
Second, we calculate the Jacobi determinant, evaluating the diagonal blocks. 

To start let us match the independent degrees of freedom present in the matrices C, D on the 
one hand and T, Q and V, U on the other hand. The matrices C and D are of size {N -\- v) x N 
and x (A^ -|- z^), respectively. Both are quaternion real, with each matrix element containing 4 real 
parameters. This leads to the following total number of degrees of freedom (d.o.f.) 

d.o±{C,D} = 8N^ + 8Ni^. (C.2) 
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On the other hand the matrix T in the parametrisation eq. ()(].H1 is triangular of size {N + v) x N: 
only Tn<i 7^ 0. The diagonal elements are of the following form |17p 



and carry only 2 real d.o.f., the eigenvalues of C. All other T„</ 7^ carry 4 real d.o.f.: 

d.o.f.jr} = 2N + 4 ^^^ ~ = 2iV^ ■ (C.4) 



The same analysis holds for the x [N + i^)-matrix Q, Qn<i 7^ 0, having N x v more non-vanishing 
elements: 

d.o.f.{Q} = 2N + 4 ^(^^~ ^) + 4^ J, = + 4iVz^ . (C.5) 
In the same way as above we define for the diagonal elements 

Qm = y* )' ^ = 1'---^' (C-6) 

containing the eigenvalues of D. The symplectic matrices U and V of size (iV + 1^)^ and N'^ a priori 
carry 2(A + z/)^ + (A^ + u) and 2A-^ + N d.o.f., respectively. However, the parametrisation eq. HC.1|) 
is not unique. Replacing the quaternion elements in all matrices by 2 x 2 blocks one can see that a 
diagonal unitary matrix U' of size 2A^ leaves the matrices T and Q triangular, keeping the eigenvalues 
unchanged: 

U UU' and V U'^V (C.7) 

where U' is U' extended to size 2(A + v) by the unity matrix. Therefore 2 A real parameters are 
redundant in eq. Furthermore, the whole lower right v x v sub-block in U can be projected 

out, subtracting a symplectic matrix with 2v'^ + v d.o.f.. We are left with V G Sp{N)/U{l)^ and 
Ue Sp{N + u)/{Sp{u)xU{l)^): 

d.o.f.{C/} + d.o.f.{y} = 2(A + uf + {N + u)-N - {2v'^ + v) + 2N'^ + N - N = 4A^ + mv (C.8) 

which together with eqs. (|C.4|) and (|C.5|) adds up to eq. HC.2|) . 

In oder to label independent d.o.f. for the computation of the Jacobian it is most convenient to 
define the following differentials: 

dC = UdCV , dC = dT + idHT - iTdJ 

dD = VUVU^ , dV = dQ + idJ Q - zQdH , (C.9) 

where 

dH = +idU^U = -iUUu 

dJ = -iVdV^ = +idVV^ (C.IO) 
Both dH and dJ are anti-self dual ^21- Their (upper N) diagonal elements are given by 



dJi,i ^ dji^i = f °, "^J' V / = 1, . . . , A . (C.ll) 



^Throughout the following we will not use summation conventions. 
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where we have used the 2N parameters of the symmetry eq. ()(].7|) to set the real parameters on the 
diagonal to zero. The remaining parameters on the upper and lower triangle are related due to the 
anti-self duality, and we only keep the latter in the following set of independent variables: 

dH : {dhi^i,dHn,i\n > I; n = 1, . . . N + i^, I = 1, . . . N} 

dJ : {dji^i, dJn,i |n > n, / = 1, . . . iV} , (C.12) 

where in dH we have also projected out the i^^-block. They label 4:N{N — l)/2 + ANu + -\-2N = 
2N'^ + ANu and 4N{N — l)/2 + 2N = 2N'^ real variables respectively, matching those of U, V in eq. 

(EIHl). 

We can now arrange the independent variables {dC, dT)} and {dT, dQ, dH, dJ} in such a way that 
the Jacobian is block triangular. It is most efficient to use quaternion real matrix elements for this 
ordering. In a second step we then have to compute the determinant of the Jacobi matrix by taking 
the product of the determinants of the diagonal blocks of quaternion elements. The variables {dC, dT>} 
are ordered row by row, starting with the lowest row of dC, until we reach its quadratic part. We then 
alternate dC and dV element by element, completing the non-quadratic part of the row with dD only. 
In detail the ordering looks as follows: 

dCN+u,l ■ ■ ■ dCN+u,NdCN+u-l,l ■ ■ ■ dCN+u-l,N ■ ■ ■ dCN+1,1 ■ ■ ■ dCN+l,N 

dCN,id'DN^i . . . dCN,Nd'DN^Nd'DN,N+i ■ ■ ■ dV^^N+u 

dCN-l,ld'DN-l,l ■ ■ ■ dCN-l^NdVM-l^NdVN-l^N+l ■ ■ ■ dVN-l,N+v ■ ■ ■ 

dCi^idVi^i . . . dCi^NdVi^^dVi^N+i ■ ■ ■ dVi^N+v ■ (C.13) 

The independent variables {dT,dQ,dH,dJ} can also be put into two full matrices without zeros, dA 
and dB of size (A^ + v') x N and N x (A^ -|- i/) respectively, and then ordered in the same fashion as 
dC and dD. The matrix dA consists in its strictly lower triangular part of dH. The diagonal part is 
composed of dhi^i and dti^i giving full quaternion real matrix elements: 

and its strictly upper triangular elements are completed by dT. Similarly we define a matrix dB by 
putting dJ in its strictly lower triangular part, the sum of dji^i and dqi^i on its diagonal: 

<i<3,, = *;). (C.15) 

and dQ on its strictly upper triangular part. These matrices dA and dB are then ordered exactly as 
in eq. ()('.13|) . or most explicitly as: 

dHfq+y l . . . dHM+u,NdHN^y^i^i . . . dHN+u~l,N ■ ■ ■ dHN+1,1 ■ ■ ■ dHN+l,N 
dHN,ldJN,l ■ ■ ■ dHN,N-ldJN,N-ldTi\f^NdQN,NdQN,N+l ■ ■ ■ dQN,N+u 

dHN-l,ldJN~l,l ■ ■ ■ dHM-l,N-2dJN~l,N-2dTN-l,N-ldQN-l,N-ldTM-l,NdQN-l,N ■ ■ ■ 
dQN-l,N+u ■ ■ ■ dH2^ldJ2,ldT2^2dQ2,2 ■ ■ ■ dT2^NdQ2,NdQ2,N+l ■ ■ ■ dQ2,N+u 

dTi^idQiA . . . dTi^NdQi,NdQi^N+i ■ ■ ■ dQi^N+u ■ (C.16) 
Next we write out explicitly the differentials in eq. ()('.9|) in terms of all independent matrix elements 
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of {dT, dQ,dH, dJ}: 

dCk,i = dTk^i + i ^ dHk^nTn,i + idhk^kTk,i — i ^ dHn,kTn 

n=l n=k+l 

l-l N 

+i ^ Tk,n dJi^n - iTk^i dji^i - i ^ T^^n dJn,i (C.17) 

n=l n=l+l 

k-1 N 

dVk^l = dQk,l + i ^ dJk,n Qn,l + idjk,k Qk,i — i ^ dJn,k Qn,l 
n=l n=k+l 
l-l N+u 

+i Qk,n dHi^n - iQk,i dhi^i - i ^ Qk,ndHn,i ■ (C.18) 



k-1 N+u 



n=l n=l+l 

For the matrix elements related through anti-self duality we have used the relation dHn.k = —dH^^ni 
where ~ denotes the conjugate quaternion ^7]. Of course it holds that dH^^n/dHf^^n is non- vanishing. 
However, we will not need to compute these matrix elements as they appear on the off-diagonal of the 
Jacobi matrix only. Let us distinguish 5 cases to show the block-triangularity of the Jacobi matrix. 
The differentiation of quaternion elements by quaternion elements used here is to be understood in a 
symbolical sense, in order to see the block structure of the Jacobi matrix. In a second step we will 
differentiate independent matrix element by independent matrix element. 

i) k>l; k = N + u,...,N + 1, l = l,...,N: 

In this case we only have differentials dCk,i simplifying to 

I 

dCk,i = i'^dHk,nTn,i ■ (C.19) 

n=l 

The variable dCk^i only depends on dH^^f 

dCk,i/dHk,i = i%i , (C.20) 

and not on variables to the right of dH^ i in the ordering eq. HC16() . that is dHk^n>u dHn<:k,rm 
or any of the dJ, dQ, dT. 

ii) k>l; k = N...,2, l = l,...,N: 

We have both differentials dC^^i an d'Dk,i simplifying to 

I N 



dCk,l — i ^ dHk^nTn,l — i ^ Tk,ndJn,l , 
n=l n=k 
I N+u 

dVk,i = i ^ dJk,n Qn,i - i ^ Qk,n dH^^i . (C.21) 

n=l n=k 

This gives the following 2x2 entry on the diagonal: 





dHk^i 


dJk,l 


dCk,i 


itl,l 


—'itk,k 


dVk^i 


-iqk,k 


iQl,l 



(C.22) 

None of the differentials depends on variables further to the right of dH^^i and dJ^^i in eq. (|C.16|) . 
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iii) k = l; k,l = l,...,N: 

k-l N 

dCk,l = dtk,k + i ^2 ^■^k,nTn,k + 'idhk^ktk,k — iik,kdjk,k — i ^ Tj^,ndJn,l , 
n=l n=k+l 
k-l N+v 

dVk^l = dqk,k + i'^dJk,nQn,l + idjk,kqk,k-iqk,kdhk,k-i ^ Qk,ndHn^i . (C.23) 

n=l n=k+l 

The diagonal block entry reads: 





{dtk,k , 


dhk,k} 


{djk,k , 


dqk,k} 


dCk,i 


1 


'itk,k 


—itk,k 





dT>k,i 





-iQk,k 


iQk,k 


1 



(C.24) 



with no dependence on variables further to the right of dT^^k and dQ^^k- 

iv) k <l] k,l = 1,...,N: 

k-l I 

dCk,i = dTk^i + i ^ dHk,nTn,i + idhk,kTk,i - i ^ dHn,kTn,i 

n=l n=k+l 

l-l N 

+i XI ^'^'^ ~ ^^'^'^ ~ ^ X ^'^'^ ^"^^'i ' (C.25) 

n=k n=i+l 

fe-1 I 

dVk,i = dQk,i + i X dJk,n Qn,l + idjk,k Qk,l - i X ^'^-^ 

n=l n=k+l 
l-l N+u 

+i X Qk,n dHi^n - iQk,l dhi^i - i ^ Qk,n dHn,i , (C.26) 

n=k n=l+l 

with little simplifications. We have 

dCk,i/dTk,i = 1 , dVk,i/dQk,i = 1 , (C.27) 
with no dependence on variables further to the right of dT]^^i and dQ^^i. 

v) k<l- k = l,...,N, 1 = N + 1,...,N + u: 

k-l N 
d'Dk,i = dQk,i + ^ X ^'^k,n Qn,l + idjk,k Qk,l — i X ^'^"■'^ ^'^-^ 

n=l n=fe+l 

/-I JV+j/ 

+i X '5'=," '^^'.n ~ "^Qk.l dhi^i - i X Ql:,n dHn,i , (C.28) 
n=A; n=l+l 

with 

dVk,i/dQk,i = 1 , (C.29) 
and no dependence on variables to the right of dQk,i- 
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This proves that the Jacobi matrix in the given orderings eqs. and (|(7.16|) is indeed block- 

triangular. It remains to compute its determinant in terms of the 1x1 and 2x2 quaternion blocks 
on the diagonal. 

In the following we collect all the contributions from the diagonal blocks. In case i) we parametrise 
for each matrix element the contributing part in eq. HC.19|) 



k,l 



r s 
t u 



dC 



kl 



a b 
c d 



rxi 

txi ux*i 



+ 



leading to 



det 



9(a, h, c, d) 



dei{(l\Bi,g{ixi,ix*i ^ixi,ix*i )) = \xi\ 



d{r, s, t, u) 

The contribution to the Jacobian from case i) for all matrix elements is thus 

N 



N N+u 

n n i-'i 

1=1 k=N+l 



n 

1=1 



Xl\ 



(C.30) 



(C.31) 



(C.32) 



In case ii) we parametrise in addition dV^^i and dH^^i as above with primed variables. We thus obtain 
for the contributing terms in eq. (|C.21() 



dCk,i 
dVkj 



a 


^) 


■ ( I'Xl 


SX*i 


c 




UX*i 


a' 








d 


d' , 


)-'{Z 





XkT Xk-S 

xp' xlu' 
Vkr VkS 



+ ... , 



+ ... . 



(C.33) 



This leads to the following 8x8 block where we only display non-vanishing elements 



det 



d{a, b, c, d, a', b' ,d , d') 
d{r, s, t, u, r' , s', t', u') 



det 



ixi 


ix* 

ixi 

ix* 


-iyk 


-iyk 


* 

-^yk 


-^yk 


-ixk 


-ixk 


iyi 


* 


iyi 


iy*i 



l^i ~ ^k\ \^l 



^2:2 



(C.34) 



Here we have used x^ilk = —zi and the following formula for determinants of matrices A' , B' ,C' , D': 



det 



A' B' 
C D' 
given by 



det(A') det(i:»' - CA'-^B'). The contribution to the Jacobian from case ii) is then 



N 



^2:2 



(C.35) 



k>l=l 



In case iii) we obtain from eq. HC.23|) for the contributing terms 
dCki = 



a b 

c d 

a' y 

d d! 



dxk idskxl - ixkdpk 
idslxk - ixldpl dxl 

dyk idpkvl ~ iykdsk 
idplvk - iyldsl dyl 



+ 



+ 



(C.36) 
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leading to the matrix 



det 



d{a, d, b, c, a' , d' , h' , d) 



d{dxk, dxl,dsk, dsl, dyk, dyl, dpk,dpl) 



det 



-1 

1 






1 








-Wk 




ixk 




* 

-^Vk 




1 






1 




-ixk 


^Vk 




-ix*k 




Wk 



1*2 2|2 
\^k ~ ^k\ 



(C.37) 



This result is obtained similarly to the previous case, with the contribution to the Jacobian now 
reading 



N 



(C.38) 



k=l 



It is easy to see that the remaining cases iv) and v) simply give a contribution unity. Collecting eqs. 
()('.32|1 . ()('.35|1 and ()('.38|1 we obtain as a final result the conjectured Jacobian in eq. 
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